Suppression of Superconductivity by Twin Boundaries in FeSe 
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Low-temperature scanning tunneling microscopy and spectroscopy are employed to investigate 
twin boundaries in stoichiometric FeSe films grown by molecular beam epitaxy. Twin boundaries 
can be unambiguously identified by imaging the 90° change in the orientation of local electronic 
dimers from Fe site impurities on either side. Twin boundaries run at approximately 45° to the 
Fe-Fe bond directions, and noticeably suppress the superconducting gap, in contrast with the 
recent experimental and theoretical findings in other iron pnictides. Furthermore, vortices appear 
to accumulate on twin boundaries, consistent with the degraded superconductivity there. The 
variation in superconductivity is likely caused by the increased Se height in the vicinity of twin 
boundaries, providing the first local evidence for the importance of this height to the mechanism of 
superconductivity. 

PACS numbers: 74.70.Xa, 07.79.Fc,61.72.Mm, 74.25.Wx 



The response of superconductivity to crystal defects 
is crucial to two forefront technological issues, namely 
the sharpness of the superconducting transition and 
the critical current. Many early studies have revealed 
a slight enhancement in the superconducting critical 
temperature Tc near twin boundaries (TBs) of certain 
conventional superconductors such as In, Sn, and Nb [1]. 
Meanwhile, TBs tend to pin vortices and so enhance the 
critical currents in the cuprate high-Tc superconductor 
YBasCuaOr-^ (YBCO) [2, 3]. The general interplay of 
TBs and superconducting properties remains unresolved. 

In the recently discovered iron-based compounds, 
the tetragonal-to-orthorhombic distortion above Tc 
typically generates a maze of TBs upon cooling 
[4], which serves as a test bed for twinning-plane 
superconductivity. Local susceptometry measurements 
with a scanned superconducting quantum interference 
device (SQUID) show an enhanced superfluid density 
along TBs in underdoped Ba(Fei_a^Coa^)2As2 [5, 6], 
compatible with SQUID magnetometry images where 
vortices avoid pinning on TBs [7]. On the other 
hand, doping-dependent TB imaging with polarized 
light, combined with bulk critical current determination 
in Ba(Fei_a^Coa^)2As2, shows a tremendous enhancement 
of critical current at the doping level where TBs are 
densest, leading to a claim that vortices are pinned 
on TBs [8]. However, in the absence of direct vortex 
imaging, this latter observation is also consistent with 
the possibility that the critical current is enhanced by 
vortex trapping between the TBs. Bitter decoration in 
Ba(Fei_a;Niaj)2As2 shows vortices clustered along lines 
in some regions of the sample, leading to a claim of 
TB pinning [9]. But the absence of simultaneous twin 
boundary imaging again leaves open the possibility that 



the aligned vortices are pinned on domains between 
parallel TBs. 

These experiments on iron-based superconductors have 
used magnetic imaging techniques, whose resolution is 
limited to approximately the penetration depth A ~ 
325 nm [10]. Because pinning may occur on the 
vortex core length scale, ^ ~ 3 nm [11], it can be 
challenging in some cases to determine from magnetic 
imaging alone whether a vortex is pinned on or near 
the TB. Scanning tunneling microscopy (STM) and 
spectroscopy (STS), which can image both TBs and 
vortices on the ^ length scale, can address this issue 
with a resolution xlOO better [12]. Additionally, in 
most cuprates and iron pnictides, chemical doping plays 
an essential role in superconductivity; thus its possible 
variation across TBs may complicate the understanding 
of the twinning-plane superconductivity [13]. As an 
alternative, the stoichiometric and structurally simple 
PbO-type /3-FeSe superconductor provides a unique 
system for addressing the variations in superconductivity 
near TBs [14]. 

Here we report on STM and STS studies of TBs in 
stoichiometric and superconducting FeSe films grown by 
molecular beam epitaxy (MBE). This allows for a direct 
probe of the superconducting order parameter near TBs 
at the nanometer length scale. Sparse Se atoms at 
near-surface Fe sites produce local dimerlike scattering 
signatures. TBs are identified by the 90° rotations of the 
electronic dimers on either side, and are seen to roughly 
orient along the diagonals of the Fe unit cells. We observe 
that (i) TBs considerably suppress the superconducting 
gap within the coherence length, and (ii) vortices tend to 
be pinned on TBs. Both observations demonstrate that 
TBs locally weaken the superconductivity in FeSe. 



All STM and STS tunneling experiments presented 
here were carried out at 4.5 K on a commercial ultrahigh 
vacuum low temperature STM apparatus (Unisoku), 
which is connected to a MBE system for in situ sample 
preparation. The base pressure for both systems is better 
than 10"^° Torr. The MBE growth of stoichiometric 
FeSe films has been described in detail elsewhere [15, 
16], and in the supplemental Material [17]. Prior to 
data collection, a polycrystalline Ptir tip was cleaned 
by electron-beam heating in ultrahigh vacuum, and 
then calibrated on a MBE-grown Ag film on a Si(lll) 
substrate. Spectroscopic measurements were made by 
disrupting the feedback circuit, sweeping the sample 
voltage, and extracting the tunneling conductance dl/dV 
using a standard lock-in technique with a small bias 
modulation of 0.1 mV at 987.5 Hz. 




bright ends positioned on two adjacent Se atoms in the 
topmost layer, suggesting that the excess Se substitutes 
into the uppermost Fe layer. Two inequivalent Fe 
positions, denoted by ji and v in Fig. 1(d), lead to 
the two orthogonal atomic dumbbells observed. Second, 
at the much larger length scale of r^ IGape-Fe — 
4.4 nm [Fig. 1(c)], the C4 symmetry is broken by 
unidirectional depressions in the density of states which 
straddle each excess Se (yellow dashes). In contrast 
to the persistence of atomic dumbbells up to 2.5 eV 
imaging bias, the larger unidirectional features exist only 
in a narrow energy range (approximately ±20 meV), 
which supports a purely electronic origin. Note that in 
Fig. 1(a) a faint stripe occurs along the upper left to 
lower right diagonal. Across this stripe, the Se-induced 
unidirectional nanostructures (electronic dimers) are 
found to rotate by 90°. This closely resembles the 
TB-induced 90° rotation of ~ Sape-Fe unidirectional 
nanostructures in slightly Co-doped Ca(Fei_a;Coaj)2As2 
[18]. We therefore argue that the observed faint 
stripe along the diagonal of Fig. 1(a) represents a TB, 
across which the a and h crystalline axes interchange. 
Here a and h correspond to the two Fe-Fe bonding 
directions, as defined in Fig. 1(e). Note that these larger 
electronic dimers always respect the crystalline a axis, 
irrespective of the stochastic distribution of atomic ji and 
V dumbbells. These observations not only provide a way 
to distinguish TBs, but also support the fundamental 
role that electronic dimers play in scattering mechanisms 
[19, 20], and by implication, transport anisotropy in iron 
pnictides [21]. 



FIG. 1. (color online) (a) STM topography of FeSe film with 
extra Se atoms appearing as bright atomic-scale dumbbells 
(y = 10 mV, / = 0.1 nA, 100 nm x 70 nm). (b) Zoom-in on 
two orthogonally oriented atomic dumbbells, labeled ji and 
h> {V — Q mV, / = 0.1 nA, 2 nm x 2 nm). (c) Larger 
zoom-in of a single excess Se atom [V — 10 mV, / = 0.1 
nA, 6 nm x 8 nm). The blue dots mark the subsurface Fe 
atoms. The depressions straddling each excess Se (marked by 
dashed yellow lines) and the TB likely stem from quasiparticle 
scattering, (d) Schematic crystal structure of /3-FeSe showing 
the inequivalent /x and v Fe sites, and (e) diagram illustrating 
a TB with Fe (blue) and Se (yellow) spheres. 

Figure 1(a) depicts a constant-current topographic 
image of an as-grown FeSe films. The localized defects 
(< 0.05%) correspond to individual excess Se atoms 
[16], which are intentionally introduced and act as 
scatterers for electrons and give rise to unidirectional 
electronic nanostructures in FeSe [15]. A more detailed 
examination shows that each excess Se explicitly breaks 
fourfold iC^) rotational symmetry in two independent 
ways, at two different length scales. First, at the atomic 
length scale, we observe two orthogonal dumbbell-like 
features, labeled as ji and v. Both atomic dumbbells are 
centered at subsurface Fe atoms [Fig. 1(b)], with their 
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FIG. 2. (color online) (a) STM topography with a TB 
indicated by a white dashed line ( F = 10 mV, / = 0.1 nA, 50 
nm X 50 nm). (b) Atomically resolved topography of a TB 
(y = 10 mV, / = 43 pA, 12 nm x 12 nm). (c) Normalized 
dl/dV spectra taken at equal separations (1 nm) along the 
white solid line (normal to the TB) in (a), (d, e) Differential 
conductance maps recorded simultaneously with image (a) at 
energies of (d) zero and (e) 2.2 meV, respectively. Tunneling 
gap is set at V — 10 mV and / = 0.1 nA. 

To find the effect of TBs on superconductivity, we use 



STS to map the superconducting gap in the vicinity of 
another TB [Fig. 2(a)]. Using the atomicahy-resolved 
STM image in Fig. 2(b), we note that the TB runs 
nearly along one of the Se-Se nearest- neighbor directions 
in the topmost layer, or equivalent ly one diagonal of the 
undistorted single-Fe unit cells. Figure 2(c) shows a series 
of differential conductance dl/dV spectra, normalized 
to the normal-state conductance spectrum above Tc 
(10 K), taken along a trajectory approaching the twin 
boundary. All curves exhibit superconducting gaps with 
clear coherence peaks. However the gap magnitude 
A, half of the energy between the coherence peaks, 
decreases when approaching the TB, suggesting that TBs 
tend to weaken the superconductivity in FeSe. This 
is further supported by dl/dV maps at zero energy 
[Fig. 2(d)] and at one of the coherence peaks at ~ 2.2 
meV [Fig. 2(e)] on the same region as Fig. 2(a). The 
TB enhances the zero-bias conductance (ZBC, inversely 
correlated with the superfluid density) and suppresses the 
coherence peaks. Our observations consistently support 
the suppression of superconductivity by TBs in FeSe. 
This contrasts with the enhanced superfluid density 
along TBs in underdoped Ba(Fei_a^Coaj)2As2 by SQUID 
measurements as well as the recent theoretical prediction 
[5, 22]. 
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FIG. 3. (color online) Superconducting gap A (black 
diamonds) and ZBC (blue circles) plotted as a function 
of distance from the TB. The blue solid line depicts an 
exponential decay, while the black one is a guide to the eye. 

Figure 3 presents the extracted superconducting gap 
A and ZBC from Fig. 2(c) as a function of distance 
d off the TB. As compared to Aq = 2.2 meV on 
TB-free regions, the superconducting gap shrinks by 
- 25% to Atb = 1.66 meV on TBs. Also, ZBC{d) 
decays with distance d from the TB as ZBC(c?) = 
ZBC(oo) + Aexp{-d/^). Here ZBC(oo) and ^ are 
the constant background and superconducting coherence 
length, respectively. Based on the exponential fitting, 
we extract a coherence length of^ = 5.5 ±0.3 nm at 
4.5 K. The coherence length ^(0) ~ 5.1 nm at zero 
temperature can be calculated from the self-consistent 
BCS gap function and ^(T) ex 1/A(T) [23] with T^ - 9.3 



K [16]. We note that anisotropic vortices have been 
recently demonstrated in FeSe and can be intuitively 
understood by direction-dependent changes in ^ with 
extrema along the a and b directions [15]. In this work, ^ 
is measured along one diagonal of the undistorted Fe unit 
cells, and thus (, ^ 5.1 nm roughly represents an average 
of £^a and ^5, which is comparable to the estimated 
coherence length of 4.5 nm from transport measurements 
[14] . Our STM images therefore demonstrate directly for 
the first time the coherence-length-scale effect of a TB on 
superconductivity in the new Fe-based superconductors. 
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FIG. 4. (color online) (a) A 150 nm x 150 nm topographic 
image with three TBs (y = 10 mV, / = 0.1 nA). 
(b) Simultaneous ZBC map showing the vortices at 2 T. 
Tunneling gap is set at F = 10 mV and / = 0.1 nA. (c) 
Schematic illustrating both TBs (white dashes) and vortices 
(red circles). 

The suppressed superconductivity and thus reduced 
superfluid density along TBs should lead to a decrease 
in energy when vortices are positioned on TBs. To 
search for such TB flux pinning, we image vortices with 
an applied magnetic field normal to the FeSe a&-plane. 
Figure 4(a) shows a topographic image with three TBs, 
where we record the ZBC map at 2 T, illustrated in 
Fig. 4(b). Previously, a pronounced ZBC peak, which 
originates from quasiparticle bound states [24, 25], has 
been found in the vortex cores of FeSe [15]. Therefore 
the yellow regions with enhanced ZBC signify individual 
isolated vortices. The observed average flux per vortex is 
~ 2.05 X 10~^^ Wb [17], consistent with single magnetic 
flux quantum, <I>o = 2.07 x 10"^^ Wb. The schematic 
depiction of vortices and twin boundaries in Fig. 4(c) 
illustrates that vortices are preferentially pinned on TBs 
[17] as long as the distance separating the neighboring 
TBs is not too large. This observation confirms that TBs 
locally suppress the superconductivity in FeSe. 

Now we consider possible explanations for the 
suppressed superconductivity by TBs in FeSe. A 
variation in chemical doping across TBs can be 
reasonably excluded [13] because the superconductivity 
develops in FeSe without any external doping, in sharp 
contrast to iron pnictides [26]. We thus consider that 
the phenomenon likely stems from the structural changes 
around TBs. Indeed, in iron-based superconductors, the 
tetrahedral geometry, both the tetrahedral angle a, and 
the anion height /ianion(pnictogen or chalcogen) above 



the Fe layer, appear to be key parameters controlling 
the superconducting transition temperature Tc [26-30]. 
For each FeSe4 tetrahedron spanning a TB, two out of 
four Se anions must be mirror symmetric with respect 
to the twinning plane [Fig. 1(e)], which will distort the 
FeSe4 tetrahedra and thus change a. However, some 
previous studies have demonstrated that a does not 
significantly affect Tc in iron chalocogenides [27, 28]. We 
therefore suggest that the tetrahedral distortion cannot 
bear sole responsibility for the observed suppression of 
superconductivity around TBs in FeSe. Then we examine 
the Se height /ise around TBs. High-pressure electrical 
resistivity measurements revealed an enhanced Tc as /ise 
is reduced [29]. In all our topographic images, up to 1 eV 
[17], TBs appear brighter than surrounding areas. This 
consistency over a wide energy range strongly suggests 
a local increase in /ise, although an electronic effect 
leading to the false appearance of increased height due 
to the STM normalization artifact cannot be completely 
ruled out. We therefore suggest that superconductivity is 
suppressed and perhaps even quenched by the increased 
hse at the TB. The well-identified superconducting gaps 
near TBs [Fig. 2(c)] may arise from the proximity effect 
between on- and off-TBs regions. 

Finally we tentatively explain the contrasting roles of 
TBs in Ba(Fei_a^Coaj)2As2 and FeSe. As has been noted 
previously, the Tc of the iron-based superconductors 
appears to reach a maximum at /lanion — 1.38 A 
[26, 29, 30]. Away from this value, Tc will abruptly 
decrease. In FeSe, h^e !::^ 1.45 A > 1.38 A [29], so the 
increased /ise must suppress Tc at TBs. However, in 
Ba(Fei_a^Coaj)2As2, /i-As — 1-34 A appears smaller than 
1.38 A [31]. Assuming that hxs increases around TBs 
of Ba(Fei_ajCoaj)2As2 as well, one can expect enhanced 
superfluid density there, in line with SQUID experiments 
[5, 6]. Moreover, /lanion may play a more important role 
in FeSe than in Ba(Fei_a^Coa^)2As2. The increased /ise 
at FeSe TBs favors the double-stripe (tt, 0) magnetic 
order, and suppresses the (tt, tt) spin fluctuations which 
are necessary for superconductivity [32]. The present 
study therefore provides evidence linking /lanion to the 
local A and thus to the mechanism of superconductivity 
in iron-based compounds. 

Our detailed STM and STS study of TBs in 
MBE-grown FeSe films has provided fundamental new 
information about the nature of superconductivity in 
iron-based materials. First, we have explicitly shown 
by direct imaging that each Fe-site impurity produces 
a local electronic dimer of size ~ 16aFe-Fe7 oriented 
along the orthorhombic a axis. Scattering from these 
dimers, although never previously directly visualized in 
real space, has been controversially suggested as the 
root cause of the transport anisotropy in iron-based 
superconductors [19, 20]. Second, we have shown by 
spatially resolved spectroscopy that TBs suppress the 
superconductivity within a superconducting coherence 



length ^. This provides a quantitative measure of the 
coherence length, ^ ~ 5.1 nm. Third, we show that 
magnetic vortices are preferentially pinned to the TBs. 
This supports the suppression of superconductivity at 
the TBs, and can inform engineering work to optimize 
vortex pinning for increased critical current. Finally, 
we show increased /ise at the FeSe TBs. This suggests 
an explanation for the contrast between TB behavior 
in FeSe vs. Ba(Fei_a^Coaj)2As2, and indeed provides 
the first local evidence for the importance of chalcogen 
or pnictogen height /lanion to the very nature of the 
superconducting mechanism in iron-based materials. 
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at Harvard University. Topographic images were partly 
processed with the WSxM software [33]. 
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We prepare high-quality FeSe films on a graphitized 
SiC(OOOl) substrate [1]. The substrate, with a resistivity of 
^ 0.1 Q-cm, consists mainly of double- layer graphene, which is 
prepared using a well-established recipe [2] . We co-evaporate 
high purity Fe (99.995%) and Se (99.999%) sources from 
standard Knudsen cells onto the substrate at 450° C. During 
the film growth, the nominal Se/Fe beam flux ratio is 20, and 
the pressure is 10~^ Torr, due to the volatile Se molecular 
beam. The high Se/Fe flux ratio compensates the losses of 
volatile Se molecules, and leads to stoichiometric FeSe films. 
The excess Se will desorb and cannot be incorporated into the 
stoichiometric FeSe due to the high substrate temperature. 

The Fei-ccSei+cc samples (where x refers to the doping 
level in the top Se-Fe2-Se layer) were prepared by depositing 
additional Se onto the stoichiometric FeSe films held at 
the lower substrate temperature of 220° C, where the excess 
Se were found to incorporate into the topmost Fe sites. 
Precise control of Se dosage enabled us to grow Fei-^Sei+x 
films at various x up to x = 10%. We found that the 
superconductivity was completely destroyed at a doping level 
of approximately x — 2.5%. In the present study, the 
small X < 0.05% barely disturbs the superconductivity. For 
example. Fig. 1(a) contains 21 excess Se atoms in a 70 nm x 
100 nm field of view, corresponding to x = 0.02%; Fig. 2(a) 
contains 5 excess Se atoms in a 50 nm x 50 nm field of view, 
corresponding to x = 0.014%; Fig. 4(a) contains 99 excess Se 
atoms in a 150 nm x 150 nm field of view, corresponding to 
X = 0.032%. 

Figure 1 shows the vortex density in various magnetic fields 
(1, 2, 4 and 8 Tesla). The density depends linearly on the field 
and agrees well with the theoretical expectation. 
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FIG. 2. (a) STM topography of a TB ( F = 10 mV, / = 0.1 
nA, 40 nm x 40 nm), and (b) the simultaneously acquired 
ZBC map showing the vortices at 2 T. The tunneling gap is 
set at y = 10 mV and / = 0.1 nA. The blue and black dashes 
indicate the TB and the two dark trenches straddling the TB, 
respectively. 



Figure 3 shows the STM topography of a TB at high voltage 
(IV). Even at high voltage, the TB appears brighter than 
the surrounding plane, indicating a true geometric height 
elevation rather than a purely electronic effect. 
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FIG. 3. (a) STM topography (25 nm x 25 nm) of a TB (blue 
arrows) at high voltage ( 1/ = 1.0 V, / = 48 pA). (b) Average 
linecut across the TB. 
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FIG. 1. Magnetic field dependence of vortex density in FeSe. 



To justify the exact position of vortex cores, we 
simultaneously imaged the high-resolution vortices and STM 
topography, as shown in Fig. 2. Here, it is clear that the 
vortex cores are located in the bright central region, rather 
than the two dark trenches. 
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